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Abstract 

Discussed is a model of the two-dimensional amnely-rigid body with 
the double dynamical isotropy. We investigate the systems with poten- 
tial energies for which the variables can be separated. The special stress 
is laid on the model of the harmonic oscillator potential and certain an- 
harmonic alternatives. Some explicit solutions are found on the classical, 
quasiclassical (Bohr-Sommerfeld) and quantum levels. 

1 Introduction 

The mechanics of an affinely-rigid body was discussed in various aspects in 
[2]> [H]-[H]- ^ n t ms paper we intend to investigate qualitatively the doubly- 
isotropic dynamical models in two dimensions, having in view applications in 
macroscopic elasticity and the theory of molecular vibrations. We expect also 
applications in dynamics of nanotubes; more precisely, we mean vibrations of 
their transversal cross-sections. On the classical level our models are completely 
integrable and may show some degeneracy properties following from hidden 
symmetries. In the two-dimensional theory there exists a relatively wide class of 
isotropic potentials which admit analytical calculations based on the separation 
of variables method [T|>J US] . In this paper the special stress is laid on the model 
of the harmonic oscillator potential and certain anharmonic models. The action- 
angle analysis and discussion of degeneracy as well as the quasiclassical Bohr- 
Sommerfeld quantization are also presented. Next we discuss the Schrodinger 
quantization procedure for such an object. We follow the standard procedure 
of quantization in Riemannian manifolds [5], i.e. we use the L 2 -Hilbert space of 
wave functions in the sense of the usual Riemannian measure (volume element). 
Some explicit solutions are found using the Sommerfeld polynomial method 
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Certain of our ideas are somehow related to those underlying the papers 

mm- 

2 Geometric description of the affinely-rigid body 

We are given two Euclidean spaces (TV, U, 77) and (M,V,g), respectively the 
material and physical spaces. Here N and M are the basic point spaces, U 
and V are their linear translation spaces, and ij £ U* ® U*, g € V* ® V* are 
their metric tensors. The space N is used for labelling the material points, and 
elements of M are geometric spatial points. 

The configuration space of the affinely-rigid body 

Q := Afl(JV, M) 

consists of affine isomorphisms of N onto M. The material labels a £ N are 
parametrized by Cartesian coordinates a K (Lagrange variables). Cartesian co- 
ordinates in M will be denoted by y l and the corresponding geometric points by 
y. The configuration $ € Q is to be understood in such a way that the material 
point a e N occupies the spatial position y — $(a). 

Let ~p denote the co-moving (Lagrangian) mass distribution in N; obviously, 
it is constant in time. Lagrange coordinates a K in N will be always chosen in 
such a way that their origin a = coincides with the centre of mass C: 

J a K dp(a) = 0. 

The configuration space may be identified then with M x Ll(U, V), 

Q = ML(N, M) ~ M x U(U, V) = M x Q int , 

where LI(f, V) denotes the manifold of all linear isomorphisms of U onto V. 
The Cartesian product factors refer respectively to the translational motion (M ) 
and the internal or relative motion (Ll(U,V)). The motion is described as a 
continuum of instantaneous configurations: 

$(i,a) 4 =cf> i K (t)a K + x i (t), (f) 

where x(t) is the centre of mass position and <p(t) tells us how constituents of 
the body are placed with respect to the centre of mass. The quantities (x\ <j) % k) 
are our generalized coordinates. 

Obviously, if we put U = V = R n , then Q int reduces to GL(n,M) and Q 
becomes the semi-direct product M. n x s GL(n,R); R" is then interpreted as an 
Abelian group with addition of vectors as a group operation. 

Inertia of affinely-constrained systems of material points is described by two 
constant quantities: 

m= dp(a), J KL = / a K a L djl(a), 
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i.e. the total mass m and the co-moving second-order moment J € U ®U . More 
precisely, it is so in the usual theory based on the d'Alembert principle, when 
the kinetic energy is obtained by summation (integration) of usual (based on 
the metric g) kinetic energies of constituents |16]-|21j. 

Substituting to this general formula the above afhne constraints Q we obtain: 

rp _ rp , rp m dx i dx j 1 d4>\ dxf B AB 
* ~~ ~2~ ~dt ~dt 2 ~dt dT • 

Certainly, if we analytically identify U and V with R™ and L1(U, V) with 
GL(n,R), then 

Tint = ^Tr (j> T ^J 



3 Some two-dimensional problems 

Now, let us discuss the two-dimensional affinely-rigid body. Considered is a dis- 
crete or continuous system of material points subject to constraints according 
to which during any admissible motion all affine relations between constituents 
of the body are invariant (the material straight lines remain straight lines, their 
parallelism is conserved, and all mutual ratios of segments placed on the same 
straight lines are constant). The conception of the affinely-rigid body is a gen- 
eralization of the usual metrically-rigid body, in which during any admissible 
motion all distances (metric relations) between its constituents are constant 
[2 . We do not take into account the motion of the centre of mass. When 
translational motion is neglected, the configuration space Q may be analyti- 
cally identified with the linear group GL(2,M), i.e., the group of non-singular 
real 2x2 matrices. The most adequate description of degrees of freedom is that 
based on the two-polar decomposition of matrices: 

4> = ODR T , (2) 

where 

cos ip — sin ip 
sin ip cos ijj 

This decomposition is connected with the algebraic Gram-Schmid orthogonaliza- 
tion. It is also know in literature as the " singular value decomposition" . The ma- 
trices 0,R e SO(2,R) are orthogonal (O t O = R T R = Id, detO = deti? = 1), 
D is diagonal and positive. The orthogonal group SO(2,R) is a commuta- 
tive group of plane rotations. Spatial rotations are described by the action of 
SO(2,R) on GL(2,R) through the left regular translations, material rotations 
are represented by the action of the rotation subgroup through the right mul- 
tiplication. In the non-degenerate case (D\ ^ -D2), the decomposition ([2]) is 



O = 



cos if — sin ip 
sin if cos f 



D = 



D 1 
D 2 



R = 
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unique up to the permutation of the diagonal elements of D accompanied by 
the simultaneous multiplying of O and R on the right-side by the appropri- 
ate special orthogonal matrices (ones having in each row and column zeros but 
once ±1 as elements). This implies that the potential energy of doubly isotropic 
models depends only on D and is invariant with respect to the permutations 
of its nonvanishing matrix elements |16j . The deformation invariants D±, D 2 
are important mechanical quantities. They are scalar measures of deformation, 
i.e. tell us how strongly the body is deformed, but do not contain any informa- 
tion concerning the orientation of deformation in the physical or material space. 
The orthogonal matrices O and R describe the space and body orientations of 
the strain. Incidentally, let us mention that the complexification of GL(2,R) 
to GL(2,C) and then the restriction to the other, completely opposite (because 
compact), real form U(2) sheds some light on our model and establishes also 
certain kinship with the three-dimensional rigid body. 

We shall consider only highly symmetric model, where J is isotropic, i.e., 
its matrix has the form /xl, \i denoting a positive constant, and I is the 2x2 
identity matrix. The isotropic kinetic energy is as follows: 



T 




dDi 
~~dT 



~dt 



4L>iA> 



dip dip 
~dt~dt 



(3) 



The matrices O and R do not enter into this equation, hence the angles 95, ip 
are cyclic variables. In these coordinates the Hamilton- Jacobi equation is non- 
separable even in the interaction-free case. However, the separability becomes 
possible in new variables, obtained by the 7r/4-rotation in the plane of the 
deformation invariants D\, D2 and by an appropriate modification of the angular 
variables. Thus, we introduce the following new coordinates: 



D 2 



D 2 ) , T] = cp — ip, 7 = ip + ip. 



In the macroscopic, phenomenological elasticity theory Di > , D 2 > 0, thus, 
a > 0, |/3 1 < a. However, describing discrete or finite systems of material points 
(e.g. molecules), one can admit singular and mirror- reflected configurations. 
Then, to some extent D\, D 2 , a, j3 may be arbitrary. The kinetic energy 
becomes then 



T 



or 



da V 
~dt) 



~dt 



(4) 



This form is both diagonal and separable. The classical Stackel theorem leads 
to the following general form of separable potentials: 



V(<p,il>,a,P) = 



V„{<p-il>) , V~(<p + tl>) 



+ V a (a) + V fi (fi). 



(5) 



4 



In this formula V n , V~ n V a , Vp are arbitrary (but regular enough) functions of a 
single variable (indicated as an argument). We consider doubly- isotropic models 
in which the potential energy does not depend on variables ip, ip (equivalently ij, 
7), i.e. V v = and V 1 = 0. Performing the Legendre transformation we obtain 
the corresponding Hamiltonian H — H a + Hp in the form: 

where p v , p^, p a , pp are the canonical momenta conjugate to ip, ip, a, /3, 
respectively, and 

The quantities H a ,Hp,p v ,p^, form a Poisson-involutive system of constants of 
motion. 

The stationary Hamilton- Jacobi equation has the following form: 

_L _L\ /V— V (—Y\ (— 1 ^ & S 

4a 2 + 4/3 2 J ^U^/ + UVV / + \2/? 2 2a 2 ) dtpd^ { > 

+ (%P\ + (tf) = ^ ( E (y a {a) + w))) , 



da J \dp / 

where E is a fixed value of the energy. Due to the fact that the variables ip, 
have the cyclic character, we may write: 

S = S v {<p) + S^WO + S a {a) + Sp(j3) =a<p + &V + S a (a) + Sp{fi) 

and the action variables are as follows: 



/ - = f 0,*l: = 277,,. •/,, = ± f \j2fi (E a - V a (a)) - ^\J^ da ' ( 9 ) 



J4 = fpi,dil> = 2*b, Jp = ± j> J 2fi {E P - Vp {fi)) - {J ^Jp dp, (10) 



where E a , Ep, a, b are separation constants. 

Remark. Let us observe that the isotropic kinetic energy 



2 

may be simply written as 

A 1 l-i 



T=>to) (11) 



T=|(i 2 + y 2 + i 2 + ? i 2 ), (12) 
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where x,y,z,u are simply the matrix elements of <j>, 



x y 
z u 



(13) 



This is formally the expression for the material point with the mass p in R 4 or the 
quadruple of such material points in R. However, in the mechanics of deformable 
bodies these generalized coordinates are not very useful for dynamical models. 

It is both convenient and instructive to use also other generalized coordinates 
in the affine kinematics. We mean coordinates in which the problem is separable; 
as mentioned, the separability in various coordinates corresponds geometrically 
to some degeneracy of the problem. And besides, those coordinates suggest some 
modifications of the potential V leading to new models of deformative dynamics, 
more realistic than the harmonic oscillator and at the same time admitting also 
some analytical treatment. As expected, in doubly isotropic models the most 
natural candidates are to be sought among orthogonal coordinates on the plane 
of the deformation invariants (D\, D-i). The most natural of them are just the 
variables a, f3 introduced above: they are obtained from Di, D2 by the rotation 
by 7r/4 in R 2 . Together with the modified angular variables 77, 7 they provide a 
system of T-orthogonal coordinates in R 4 , i.e., in the space of variables x, y, z, 
u. To be more precise, they are orthogonal coordinates for the metric element 
dx 2 + dy 2 + dz 2 + du 2 on which the kinetic energy T is based. And moreover, 
as said above, they are the nice separation variables for T in the Stackel sense. 
Other natural T-separating variables are obtained as some byproducts of a, /3. 
The most natural of them are polar variables in the R 2 -plane of the pairs (a, (5). 
In certain problems it is analytically convenient to use the modified "polar" 
variables r, § given by 

r * r r ■ d 

a = y'rcos— 1 p — y'rsm—. 
Obviously, the "literal" polar variables p, e are defined by 

•d 

a = pcose, j3 — psine; p=^/r, e = — . 

The natural metric on the manifold of 2 x 2 matrices, 

ds 2 = Tr (d<p T d(j)) = dx 2 + dy 2 + dz 2 + du 2 , 



becomes then 



o o o o *t/ o 1 o o 

ds = rcos — drj +rsin — d-i + —dr + -dti 
2 2 4r 4 

= dp 2 + p 2 de 2 + p 2 cos 2 e drj 2 + p 2 sin 2 e dj 2 

= dp 2 + -p 2 d-d 2 + p 2 cos 2 - dr? + p 2 sin 2 - d-j 2 . 
4 2 2 



G 



Evidently, kinetic energy is then expressed as follows 



T = 




The above crowd of expressions is due to the fact that different conventions 
are better suited to different analogies: the two-dimensional homogeneously de- 
formable body and three-dimensional spherical top with dilatations. Physically 
we are interested here in the first problem, however, certain aspects of the sec- 
ond one (spherical top with dilatations) are formally useful and the mysterious 
link between them is interesting in itself. 

Let us notice that (r, §) may be interpreted as polar coordinates in the two- 
dimensional space of quantities 2D\D 2l D\ 2 — D 2 2 , 

2L>iL> 2 = r cos 1?, D x 2 - D 2 2 = r sin tf, (14) 

or, inverting these formulas, 

2 7^ 2 , r-> 2 x a i_„/r._\ l/^i D-2 



r = p z = D 1 z + D 2 \ tani? = tan(2e) = -^--^J. (15) 

Therefore, d refers to the shear degrees of freedom, whereas r — p 2 is some kind 
of the measure of size. More precisely, dilatation is measured by the product 
DxD 2l thus, 

.... 

r = — (16) 

cos V 

contains an "admixture" of the shear parameter Nevertheless, just like D\D 2 
it is a homogeneous function of degree 2 of (Di, D 2 ). The shear parameter 1? is 
evidently a homogeneous function of degree zero. 

It is also convenient to parametrize deformation invariants as follows: 

n ( a + b \ n 

D\ = cxp — - — , L> 2 = exp 



Then 



a = 4=(Di +D 2 ) = V2et cosh^, (3 = -^=(D 1 - D 2 ) = V2et sinh^, 
\/2 2 y2 2 

L>iL> 2 =e a , Di 2 + D 2 2 = 2e a cosh6, D x 2 - D 2 2 = 2e a sinh b, ^=e b , 



sin$ = tanh&, cos-d = — tani9 = sinhfc. 



D 2 

tan ?9 = sinh ft. 

cosh b 
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These simple formulas shed some light onto the link between two-dimensional 
homogeneously deformable body and three-dimensional top. Nevertheless, this 
link is still rather mysterious and obscure. 

For the completeness let us also mention about other orthogonal coordinates 
on the plane of deformation invariants: 

(i) Elliptic variables (k, A), where 

a = \/2coshKCOS A, (3 = v^sinhftsinA. 



(m) Parabolic variables where 



(Hi) Two-polar variables (e, /), where 

c sinh e c sin / 



a = 



cosh e — cos / ' cosh e — cos / ' 



and c is a constant. 

For our analysis of the deformative motion the parabolic (£,5) and two- 
polar variables (e, /) are non-useful, because the corresponding Hamilton- Jacobi 
equations are non-separable even in the non-physical geodetic models, i.e., ones 
with vanishing potentials. In the elliptic coordinates (k, A) the metric underlying 
the kinetic energy takes on the form: 

ds 2 = Tr (d(f) T dcj>) — (cosh 2 k — cos 2 A) dn 2 

+ (cosh 2 k — cos 2 A) dA 2 + cosh 2 k cos 2 Xdrj 2 + sinh 2 n sin 2 Acfy 2 . 

The general Stackel-separable Hamiltonians H = T + V in the variables 
(a, (3, rj, 7), (r, 1?, 77, 7) and (k, A, i], 7) have respectively the form: 
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+ Va (a) + V^) + ^ + ^, (17) 
// = — i Irp, 2 -f - ( ^ W + + , ,: 



a 2 /3 2 
J_ 2 + 1 / Py 2 + Py> 2 + 2p y py, cos t? 

2/J \ r r \ sin 2 1? 

+ Wr) + ^) + J^ + JiW, (18) 

r r cos 2 f r sin I 

H = 1 f p * 2 1 

4/i I (cosh 2 k — cos 2 A) (cosh 2 k — cos 2 A) 



Pn 2 Pi 2 



cosh 2 k cos 2 A sinh 2 k sin 2 A , 

V k (k) Vx(X) 



+ 
+ 



— w + — (19) 



2 (cosh 2 k — cos 2 A) 2 (cosh 2 k — cos 2 A) 
2 cosh 2 k cos 2 A 2 sinh 2 k sin 2 A 



Let us observe that, obviously, 

cosh 2 k — cos 2 A = sinh 2 k + sin 2 A 

and it does not matter what is written in the corresponding denominators above. 
Making use of this fact we immediately see that when the problem is doubly 
isotropic, i.e., V v , V 1 are constant, then obviously (p v ,p 7 ), equivalently (p v ,p<p), 
are constants of motion but also there is a separation of the Hamilton-Jacobi 
equation in the variables k, A. Therefore, there are two additional constants of 
motion and the problem is integrable. Those constants of motion are given by 

h K cos 2 A — h\ cosh 2 k h K cos 2 A — h\ cosh 2 k 



L = 



2 (cosh 2 k — cos 2 A) 2 (sinh 2 k + sin 2 A) 
h K sin 2 A — h\ sinh 2 k h K sin 2 A — h\ sinh 2 k 
2 (sinh 2 k + sin 2 A) 2 (cosh 2 n — cos 2 A) 



where the auxiliary quantities h K , h\ are not constants of motion and are re- 
spectively given by 

1 (_ 2 , ^ \{P v ~Pi>) 2 + 2fiV K \{p v +P*) 2 + ^V K 



l[i \ cosh k sinh k 

1 p\ + 2\xV\ 1 



2^ V cos 2 A sin A 

we remember that V K , V\ are constants here. 
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Therefore, we have the involutive system of constants of motion (their Pois- 
son brackets do vanish), and 

H = K + L 

has the vanishing Poisson brackets with all of them, i.e., with p v , p^ (i.e., with 
p v , p 7 ), K, L. 

The elliptic coordinates and the corresponding separable models are not very 
interesting for applications. From this point of view the "polar" coordinates 
(r, 1?), or equivalently (p, e), are much more useful. The configurational metric 
tensor is then expressed as follows: 

1 r 

ds 2 = Tr (d(/) T d<j)) = —dr 2 + -d-d 2 + rdtp 2 - 2r cos tfdipdip + rdtp 2 

= dp 2 + ^p 2 (rfi? 2 + d(2(p) 2 - 2 cos dd{2ip)d(2iP) + d{2iP) 2 ) 

= -!- (dr 2 + r 2 (dS 2 + d<P 2 -2 cos 6d$d* + d^ 2 ) ) , 

where, obviously, the doubled angles are used, © = $,$ = 2ip, ^ = 2tp. This 
expression is very interesting in itself. We used here three alternative systems of 
symbols, each of them convenient and suggestive in some areas of applications. 
It is seen that the expression 

da 2 = dQ 2 + - 2 cos 9d$d* + rf* 2 

is exactly, up to a constant multiplier, identical with the doubly- invariant (i.e., 
both left- and right-invariant) squared metric element on the rotation group in 
three dimensions, SO(3,R), or on its covering group SU(2). This identification 
is based on interpreting $, 0, ^ as Euler angles. More precisely, to be literal 
in this analogy, one should change the sign at then one obtains the usual 
expression 

da 2 = de 2 + d$ 2 + 2 cos Qd<f>d^ + d^ 2 . 
This metric underlies the kinetic energy expression for the spherical top, 

m I f fdQ\ 2 fd<S>\ 2 n ^d$d^ fd^>\ 2 \ 
T =2[U) + U) +2COS9 ^ + U) J' 

In mechanics of gyroscopic systems $, 0, ^ are referred to respectively as the 
precession, nutation and rotation angles. This, of course, has nothing to do with 
our object, i.e., homogeneously deformable two-dimensional body; such a body 
has only one rotational degree of freedom. The analogy is formal, nevertheless 
instructive and effective in the computational sense. The idea has to do with 
the "concentric" parametrization of K 4 . As mentioned, the Cartesian variables 
x, y, z, u, i.e., matrix elements of the configuration matrix <j) 1 are non-effective 
when investigating deformations. This was just the reason to use the two-polar 
decomposition and the corresponding coordinates (D\, D2, <p, ip) or (a, j3, tp, ip). 
The two "radii" (Di,D 2 ) or (a, (3) have to do with the purely scalar deformation; 



10 



(<p,ip) (equivalently (77,7)) are angular variables of compact topology (orienta- 
tion of deformations in the physical space and in the body) . The " concentric" 
parametrization consists in encoding the possibility of unbounded motion in the 
radial variable in R 4 , 

p = sfr = v/s 2 + y 2 + z 2 + u 2 = v^Di 2 + D 2 2 = Y^Tr(</) T (/>) = Vt^G, 

where the symbol G is used for the Green deformation tensor expressed in 
the Cartesian coordinates. More geometrically, we are dealing here with the 
deformation invariant: 

P = Vv ab Gab = \j gijft a4>' j bV AE 

g, 77 denotes respectively the spatial and material (reference) metric tensors. 

Degrees of freedom orthogonally transversal to the radial variable p (or equiv- 
alently r) describe the geometrically bounded aspect of motion. Those modes 
of motion are encoded in the concentric spheres in R 4 , in particular, in the unit 
sphere given by equation p = 1, i.e., r = 1. But it is well-known that the group 
SU(2), i.e., the group of unitary unimodular matrices and the covering group 
of SO(3,R), may be naturally identified with the unit sphere 5 3 (0,1) C R 4 . 
And in this way this sphere may be parametrized with the use of the Euler 
angles G, The parametrization of R 4 with the use of variables (p, $, Q, 
or (r, $, 9, Vt) is rather nonusual, however well-suited to the description of the 
three-dimensional rigid body with imposed dilatations or, as we see, to the 
description of the two-dimensional homogeneously deformable body. In other 
applications one uses rather spherical systems of coordinates in R 4 , e.g., r, A, 
/i, v, where 

x = r sin A cos /x cos u, 

2 • \ 

x = r sin A cos p, sin v, 

x 3 — r sin A sin /j, 
x A = r cos A. 

Let us mention that the isotropic harmonic oscillator may be described ob- 
viously in terms of those variables, and the expression of Hamiltonian through 
the action variables J r , J>, J M , J^, in analogy to (|32l) below, is given by 

H = tj(2J r + J x + Jf, + J„), (20) 

where the degeneracy, i.e., the resonance between J r , J\, J M , J„ is explicitly 
seen. 

One can also use certain mixed type parametrizations in K 4 , e.g., represent- 
ing it as R 3 x R, I 2 x R 2 and taking spherical coordinates in M. 3 or polar ones in 
one or two copies of R 2 . In all such coordinate systems the isotropic harmonic 
oscillator is separable and this is some aspect of its very high, total degeneracy. 

However, it is hard to realize a wider class of realistic applications of these 
coordinates, e.g., in elastic and similar problems. Unlike this, the apparently 
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exotic parametrization in terms of the "radial distance" p and "Euler angles" 
<£>, 0, ^ offers certain models of potentials which are both separable and quali- 
tatively physical. 

We have quoted the general Stackel-separable Hamiltonian in the variables 
(r, §, <p, -0) {H]). It is doubly isotropic when the shape functions V v , V 1 are put 
as constants. Obviously, the corresponding terms V n j cos 2 ($/2), V 7 / sin 2 ($/2) 
may be simply included into Vg(-d). We have the following four constants of 
motion in involution, responsible for separability: 

• Pip,Prp, i-e-, equivalently p,,,p 7 , 

• h ® = -k^S (Pf 2 +P4> 2 +2p^cosi?) + |p tf 2 + V#(tf), 

• H = T + V = H r + ^-, where, however, the two indicated terms in if, 
namely 

H r = -rp r 2 + V r (r), ^ 
fi r 

are not constants of motion when taken separately. 

The term V r stabilizes the radial mode of motion which without this term would 
be unbounded, therefore physically non-applicable in elastic problems. The 
term V$ is responsible for the shear dynamics. Let us stress that in spite of the 
"angular" character of § the shear mode of motion is also non-compact. It is 
just seen from the fact that the shear is algebraically expressed by the quantity 
tan??, which is unbounded. Therefore, in certain problems some non-constant 
expression for V& is also desirable. Even if we use V r proportional to r = p 2 , 
any model with non- vanishing Vg introduces some anharmonicity. Particularly 
interesting is the following simple model: 

V = V ri r ) + YM = C r+ M =C (-^ + D i 2 + D > 2 ). (21) 
W r 2 rcostf \DiD 2 2 J V ' 

The model is perhaps phenomenological and academic, however, from the "elas- 
tic" point of view it has very physical properties: it prevents the collapse to the 
point or straight-line, because the term \jD\D-i is singularly repulsive there, 
and at the same time it prevents the unlimited expansion, because the harmonic 
oscillatory term G(D\ 2 + D 2 2 )/2 = C(a 2 + /3 2 )/2 grows infinitely then. There 
is a stable continuum of relative equilibria at the non-deformed configurations 
when D\ = D 2 = 1 . Expansion along some axis results in contraction along the 
perpendicular axis, because 



dD x dD 2 

at Di = D 2 = 1. This qualitatively physical potential of nonlinear hyperelas- 
tic vibrations is separable, therefore, at the same time it is also in principle 
analytically treatable. Its structure seems to suggest some three-dimensional 
models with the attractive harmonic term proportional to (Z?i 2 + D 2 2 + -D3 2 ) 
and some collapse-preventing term, e.g., one proportional to {DiD 2 D^)^ p or 
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(D 1 D 2 )~ P + (D 3 Di)~ p + (D 2 D 3 )~ P , p > 0, however, there is no chance then for 
separability and integrability. 

In the chapter below we begin with some problems concerning the harmonic 
oscillator, 

V{a,P) = £ («* + (A 2 + A. 2 ) 

= ^+9 2 | Z 2 + « 2 )^Tr(ft),C>0. (22) 
and then discuss some natural anharmonic modifications. 



4 Harmonic oscillator and certain anharmonic 
alternatives 

The expressions J Q , Jp depend on potentials V a (a), Vp(/3), respectively. After 
specifying the form of these potentials we can obtain the Hamilton function H 
as some function of our action variables, i.e., H = E(J a , Jp, J Vl J^). We can 
find the explicit dependence of the energy E on the action variables and the 
possible further degeneracy. We will also perform the usual Bohr-Sommerfeld 
quantization procedure for our model. 

Hence, we consider the model of the harmonic oscillator potential (|22|) . Some 
physical comments are necessary here. Namely, the potential (|22|) describes only 
the attractive forces which prevent the unlimited expansion of the body. Its non- 
physical feature is that it does not prevent the collapse, i.e., the contraction to 
the null-dimensional singularity. It attracts to the configuration D\ = D 2 — 
instead than to the non-deformed state D\ = D 2 = 1. Nevertheless, the 
model may be useful in some range of initial conditions. Except the subset of 
measure zero in the manifold of those conditions, the collapse to D\D 2 = is 
prevented by the centrifugal repulsion. And the collapse missbehaviour of 122|) 
is not very malicious when the system is discrete. Obviously, (fTT|) and (|22|) 
describe the isotropic harmonic oscillator in R 4 or the quadruple of identical 
one-dimensional oscillators in K. In this sense the solution is obvious and a 
priori known. Nevertheless, the model is a useful step towards investigating 
more realistic ones. And another point is very important. Namely, the very 
strong degeneracy of this model has to do, as usually, with the separability of 
the Hamilton- Jacobi equation in several coordinate systems. 

After some calculations we obtain the dependence of the energy E = E a +Ep 
on the action variables as follows: 

E = [4J + \J V - J v ,| + \J V + J |] , J = J a + J p , (23) 

47T 

where u> = y/C/n and 

E a — — (4Jo, + | J v — J^|) , 
E = ^(iJp + \J v + J^\). 
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Then performing the Bohr-Sommerfeld quantization procedure, i.e. suppos- 
ing that J = nh, J v = mh, Jw, = Ih, where h is the Planck constant and 
n = 0, 1, . . . ; m, I — 0, ±1, . . ., we obtain the energy spectrum in the following 
form: 

E = -hw[An+\m-l\ + \m + l\\. (24) 
We may rewrite this formula as follows: 
(i) if \m\ > \l\, then m 2 > I 2 and 

E = hw (2n ± m) , (25) 

(ii) if |m| < |Z|, then m 2 < I 2 and 

E = huj{2n±l), (26) 

(Hi) if |m| = then m 2 = I 2 and 

E = Tiw (2n ±m) — Huj (2n ± Z) . (27) 

And similarly, on the purely classical level of the action variables we have 
the following formulas: 

(i) in the phase space region where \ J V \ > \J^>\: 

E=^(2J±J ip ) = ^- (2J a + 2J p ± J v ) , (28) 

Z7T Z7T 

(ii) in the region where \J V \ < \J$\: 

E=^-{2J± J 4 ,) = (2J a + 2J p ± J^) , (29) 

(Hi) on the submanifold where J v = J^: 

E= ^-(2J± J v ) = ^-(2J± J+). (30) 

Z7T Z7T 

The total degeneracy of the doubly invariant model with the potential (|22l) 
is a priori obvious because in coordinates (x, y, z, u) it is explicitly seen that 
we deal with four-dimensional isotropic harmonic oscillator (equivalently-with 
the quadruple of identical non- interacting oscillators). If we use coordinates 
(Di,D 2 ,ip,ip), or equivalently (a, j3,ip, tp), then the total degeneracy is visual- 
ized by the fact that the action variables J a , Jp, J v , enter (|2"8"|) with integer 
coefficients, with the vanishing one. Similarly in (|29p they are also combined 
with integer coefficients, but now the coefficient at J v does vanish. The third 
case (|30|) is, so-to-speak, the seven-dimensional " separatrice" submanifold. The 
existence of those regions with various expressions for the functional dependence 
of energy on the action variables is due to the fact that the coordinate system 
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(Di,D 2 , <p, ijj) is not regular in the global sense and has some very peculiar sin- 
gularities. Nevertheless, it is just those coordinates that are more natural and 
physically lucid in dynamical problems. 

The quasiclassical degeneracy of the Bohr-Sommerfeld energy levels is due to 
the fact that the quantum numbers may be combined in a single one, although 
in slightly different ways in three possible ranges. Let us observe that in (|25p 
the quantum number I still does exist although does not explicitly occur in 
the formula for E. It runs the range |Z| < \m\ and labels quasiclassical states 
within the same energy levels. And analogously in the remaining cases (l2l)|) , 
(j2"Tj) . The action variables J v , and the corresponding quantum numbers 
m, I take symmetrically the positive and negative values, thus, as a matter of 
fact, the ambiguity of signs in the above formulas (|25"|) -(f2"7 f does not matter 
when the values of energy in stationary states arc concerned. Nevertheless, this 
ambiguity is essential for classical trajectories, namely, for different signs the 
orbits or rather their angular cycles are "swept" in opposite directions. 

Let us observe that the formulas (|2"8|) - (f5tJ)) resemble the action- angle descrip- 
tion of the two-dimensional isotropic harmonic oscillator in terms of usual polar 
coordinates (r, <p) on M. 2 . Namely, the Cartesian formula 

E = lo(J x + Jy) (31) 

is then alternatively reformulated as 

E = uj(2J r + J v ). (32) 

The ratio 2 : 1 of coefficients is due to the fact that the total angular rotation 
in the ip- variable is accompanied by the exactly two total cycles of " libration" 
in the r-variable. The analogy is neither accidental nor superficial. For the 
deformative motion the deformation invariants D±, Z?2, i-e., stretchings, are 
analogues to the radial variable r, whereas the two-polar angles ip, ip describing 
the spatial and material orientation of stretchings play a role similar to the polar 
angle tp in material point dynamics on R 2 (do not confuse-the symbol tp is used 
in two different meanings). This is just the reason for the 2 : 1 ratio in (j2"0|) and 

12BH3JD. 

Let us now review certain still isotropic, but anharmonic modifications of the 
harmonic model of affine vibrations (fTT]) and (|22|) . They are based on the use of 
variables (a, (3, ip, ip) or (p, i?, p>, ip). The corresponding potentials are given by 

V(a, P) = £ (a 2 + ±) + ^ = 2 (a 2 + /3 2 ) + ^, (33) 

2 V P J P 2 2 p z cos 2 | 

where in both formulas C denoting some positive constant. 
Using the former symbols we have 
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An important peculiarity of these models is that they have the stable equilib- 
ria in the natural configuration D\ = Di = 1, so they are viable from the elastic 
point of view. And both of them are separable ( (|33l) in the obvious additive 
sense), therefore, the corresponding Hamiltonian systems are integrable. 

One can explicitly calculate the action variables that correspond to and 
i.e., (J a , Jj3, Jtp, J^) and (J r , J#, J v , J^). They are some functions of the 
separation constants (one of them is the energy E) . Eliminating other constants 
one obtains the expression of E, or more precisely, of the Hamiltonian H, as a 
function of action variables. 

For the model (I3"3l one obtains 



E = i^ (4(J Q + J/t) + \J V + J+\ + y'«;W" + !•/,- - ■>- ) 
where, as usually, we denote 




It is seen that the collapse-preventing term C/a 2 in V a partially removes the 
degeneracy. Evidently, there is no longer resonance between tp and ip. The 
resonance between a and j3 obviously survives; their conjugate actions J a , Jp 
enter the energy formula through the rational combination J = J a + Jp and the 
corresponding frequencies are equal: 

u 

Vu=Vp, = —. 

7T 

We use here the standard formulas: 

dE dE dE __ dE 

a dJ a ' 13 dJp ' v dJ<p ' ^ dJtf, 

There are two phase-space regions given respectively by J v + > and J v + 
Jjf, < 0. In any of these regions there is a resonance between 7 = p + ijj and a, 
j3. This is seen from the formulas 

'J \n J' j-j ~~\~~ J*yi Jib Jn J y • 

In the mentioned regions we have respectively 



E = ^ (j 4 J " + AJ P ± 2 J 7 + \/^^ 2 C + J r , : 
This implies the following independent resonances: 

v a — vp = 0, v a =F 2z^ 7 = 

or, equivalcntly, 

v a — Vfj — 0, v/3 =F 2z/ 7 = 0. 
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Therefore, in any of the mentioned regions, where J 7 > or J 7 < 0, the system is 
twice degenerate and the closures of its trajectories are two-dimensional isotropic 
tori in the eight-dimensional phase space. 

Using the primary variables if, ip, we have the following expressions for v v , 

.. _ cj ( ±l t 2{J V -J^) ^ 



2 



2(J^ — J t 



the ± signs respectively in the regions where J v + J^ > or J v + < 0. Then, 
taking into account that 

dE 

u = TXV a = nvp = nv = — , 
we have the following degeneracy conditions: 

v a - vp — 0, v a =F 2v v =F = 0, 

respectively in the regions where J a + Jp > or J a + Jp < 0. Obviously, in the 
second equation, v a may be equivalently replaced by vp. 

The corresponding quasiclassical Bohr-Sommerfeld spectrum is given by 



E = huj Un +\m + l\ + yj(m- I) 2 + . (35) 

Another interesting model is (|34l) . separable in the variables (p, #), i.e., equiv- 
alently (r, ??). Then we obtain 



E = f 

47T 



^4(J r + J # ) + | J v + J^l + y / 64^7r 2 C+ (J v - J^) 2 



- (^4(2J p + J,,) + | J v + + y/Mfiir 2 C + (X ~ J^Y 

Again there is only a two-fold degeneracy and the system is not periodic. 
Trajectories are dense in two-dimensional isotropic tori. Degeneracy is described 
by the following pair of independent equations: 

v p - 1v® = 0, v$ =F 2z^> =F 2i^ = 0, 

respectively in the phase-space regions where J v + > or J v + < 0. 
Obviously, the second equation may be alternatively replaced by 

v p T T 4z> = 0. 
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The corresponding quasiclassical Bohr-Sommerfeld spectrum is given by 



where the quantum numbers n, m, I, refer respectively to the action variables 
J, Jtp, J^, and the system is twice degenerate. Quasiclassical energy levels are 
labelled by two effective quantum numbers, namely, {An + m + I) and (m — I), 
and there is also an obvious degeneracy with respect to the simultaneous change 
of signs of to and I . 

Let us mention that some anharmonic potentials independent of i9, e.g., the 
first term in (l34l) . are also of some practical utility as models of a bounded 
motion. The point is that, as seen in formula (1151) . the variable r depends 
both on the area of the body (its "two-dimensional volume") and on the shear 
parameter. Therefore, to be bounded in r implies to be bounded both in the 
"volume" and shear degrees of freedom. Due to the separability, the motion in 
(<p, i?, ip) -variables is geodetic in the sense of invariant metric tensors on SO (3, R) 
or SU(2). And this problem is mathematically isomorphic with the motion of 
the free spherically-symmetric rigid body in the three-dimensional space (purely 
rotational one, without translations in R 3 ). 

Another helpful model would be one with V^(i?) = Acosd, where A denotes 
some constant. The resulting problem is isomorphic with that of the three- 
dimensional heavy top. 

It is not excluded that some more general problems from the realm of three- 
dimensional gyroscopic dynamics, e.g., symmetric top, might be also of some 
mathematical usefulness when studying the two-dimensional affine motion. 

5 Quantized problems 

Classical dynamical models described above may be easily quantized in the sense 
of Schrodinger wave mechanics on manifolds. And those rigorously solvable on 
the classical level are so as well on the quantum level. 

Let us fix some notation. Let Q be a differential manifold of dimension n 
with the metric tensor G. The components of G with respect to some local 
coordinates q 1 , . . . , q n will be denoted by Gy and the components of the con- 
travariant inverse of G will be denoted by G lJ ; by definition, GifcG fcj = <V . The 
determinant of the matrix [G^] will be briefly denoted by the symbol |G| (no 
confusion between two its meanings); it is well-known that, this determinant is 
an analytic representation of some scalar density of weight two; the square root 
\/jG] is a scalar density of weight one. The invariant measure induced by G 
will be denoted by ju; analytically its element is given by 



Operators of the covariant differentation induced in the Levi-Civita sense by 
G will be denoted by Vj. The corresponding Laplace-Beltrami operator A is 




dfl(q) = y/\Gyj\dq 1 ---dq n . 
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analytically given by 



or explicitly, when acting on scalar fields, 

1 



A* 



y- 

p) n i 



\G\ *r! dq 

1 1 «J 



|r;|G V 



SE' denoting a twice differentiablc complex function on Q. 

Wave mechanics is formulated in L 2 (Q,J1), the space of square- integrable 
functions on Q with the scalar product meant as follows: 



(*|*) := / *(g)*(«)d£(g). 



The operator A is symmetric with respect to this product, and Vj are skew- 
symmetric. The metric G underlies the classical kinetic energy, therefore, the 
classical energy /Hamiltonian function 



becomes the operator 



H = A 

2fi 



V. 



Then, denoting and ordering our coordinates q % as (<p, ip, a, 0) in the Cartesian 
case we have for explicitly separable isotropic potentials: 



[Gi. 



a 



(3 2 /3 2 







a 



a 2 +p 2 



1 
1 



where 



H« = 



H = H a + Hn 



1 



2n 



A + V(a,/3), 



(36) 



(37) 



d 2 

da 2 
i 



ld_ 

a da 



+ V a (a), (38) 



and S = (h/i)d/dip is the spin operator, the generator of spatial rotations about 
the current spatial position of the center of mass, whereas E = (h/i)d/dip is the 
"vorticity" operator, the generator of material rotations. Operators H a , Hp, S, 
S are the quantum constants of motion. They also commute with each other 
(they represent co- measurable physical quantities). 
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Those formulas follow from the expression of A in coordinates (tp, ip, a, 0) 

9 2 * 9 2 * 1 9* 19* / 1 1 \ /9 2 * 9 2 * 
A* = 1 1 1 h 1 1 

9a 2 B0 2 aBa 80 \ia 2 A0 2 ) \ dip 2 dip 2 

1 1 \ 9 2 * 



20 2 2a 2 J Bydip 



(40) 



Separable solutions of the stationary Schrodinger equation if* = B* have the 
form: 

9((p,^a,P) = f v (<p)UW)f a (a)fp(p), (41) 

where f v {<p) — e™ v , fy(tp) — e tl ^ (m, I are integers) and f a (a), fp{0) are the 
dcformative wave functions. 

Hence, the stationary Schrodinger equation with an arbitrary potential 

V(a,0) = V a (a) + Vp(0) 

leads after the standard separation procedure to the following system of one- 
dimensional eigenequations: 

d 2 f a (a) , 1 df a (a) (m-tf 2/x 

~~ da 2 a - da 4a 2 — ' ' ft 2 ^ ►<«(«)) /«( a ) = °> ( 42 ) 

_ ^ 2_ + ^ i^-^^ + F ( ^ - W)) = °" (43) 

It is natural to expect that for potentials ([5]) the resulting Schrodinger equa- 
tions should be rigorously solvable in terms of some standard special functions. 
The most convenient way of solving them is to use the Sommerfeld polynomial 
method [13j [14]. In this method the solutions are expressed by the usual or 
confluent Ricmann P-functions. They are deeply related to the hypergeomet- 
ric functions (respectively usual F% or confluent F2). If the usual convergence 
demands are imposed, then the hypergeometric functions become polynomials 
and our solutions are expressed by elementary functions. At the same time 
the energy levels are expressed by the eigenvalues of the corresponding oper- 
ators. There exists some special class of potentials to which the Sommerfeld 
polynomial method is applicable. The restriction to solutions expressible in 
terms of Riemann P-functions is reasonable, because this class of functions is 
well investigated and many special functions used in physics may be expressed 
by them. There is also an intimate relationship between these functions and 
representations of Lie groups. 

Let us now quote some formulas for quantized problems separable in coordi- 
nates (r, 0, \P) (equivalently (p, <I>, O, *)), namely, the quantum counterparts 
of classical models (|18[) . One can easily show that the Laplace operators take 
on the form: 

9 2 * 9* 1 /9 2 * „ 9 2 * 9 2 * 

A* = Ar— T +%—+ + 2 cos ??—— + — 3- 

ar z or r sin v \ 9<^ z o<poip atp^ 

, /9 2 * 9* 
- | 9^ + COt ^ 
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i.e., 

<9 2 * 3 9* 4 /a 2 * ^ a 2 * <9 2 * 

A * = -rr^r + — ^ + ^ 9— tttt + 2 cot 9 — - T + — 7 

dp 2 p dp p 2 sin 2 9 V d ® <9* 2 

4 /a 2 * ^a* 

We assume the doubly isotropic separable potential energy (|2"Tj) . i.e., 

The corresponding Schrodinger equation separates and, taking into account the 
cyclic character of angular variables (p,ip, we put 

9(<p,i/>,r,#) = e m *e a ^f r (r)h{d) = e im * ^ f p {p) U{§) , (44) 

where m, I are integers. 

Quantum integration constants responsible for this separability are given by 
operators: 

• = 1-L = S- spin, 



A A. 

i dtp 



P'P = Ilk = V- vorticity, 



• = 2^FH (Pv 2 + 2cosi?p^ +p^ 2 ) cottf^) + Vg, 

• H = H r + H. e = H r + \hi) = H p + jjhi}- energy, 
where the " radial energy" is given by 

The four mentioned constants of motion p v , p\p, h&, H are pairwise commuting 
and therefore they represent co-measurable physical quantities. 
Warning: the two indicated contributions to H, i.e., H r and H$ = h$/r are 
not constants of motion. 

The stationary Schrodinger equation for the factorized wave function (|44l) 
reduces to the following pair of ordinary Schrodinger equations (Sturm-Lioville 
equations) for the factors depending only on one variable, respectively 1? and r 
(or p): 

h#U = e#U, (45) 

Hrfr + —fr=Ef r , i.e., H p f p + ^f p = Ef p . (46) 
r p* 

The procedure is first to solve the ^-equation and then to substitute the 
resulting eigenvalues e# to the r/p-equation. Then one obtains (at least in 
principle) the energy levels E. 
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It was mentioned that there exists some strange relationship between the 
two-polar parametrization of GL(2,R) and the Euler angles and scale parame- 
ters of rigid body with dilatations. There is some very interesting aspect of this 
link, which we noticed first quite accidentally, on the purely analytical level, 
before the trivial geometric meaning of this surprise became evident to us. This 
artificial detour (wandering about) was due to the fact that by chance we in- 
vented our separating coordinates (r, i?) better (p, i?) just where they are rather 
obscurely hidden, namely as polar parametrization of the pair of quantities 
{2D l D 2 ,D 1 2 - D 2 2 ) 

Namely, differential eigenequations (|4"5l) , (|4^|) may be explicitly written down 
as follows: 

4^ 8 f+!(M*+T))^ <«> 

where to, I are integers in ^ as coefficients at the angles ip, ip in complex 
exponential functions (eigenfunctions of p v , fL,). Let us now divide by 4 the 
nominator and denominator in the bracket expression (|47[) and formally admit 
half-integer coefficients. We can rewrite our equations as follows: 

where now the numbers to, Z are assumed to run over the set of non-negative 
integers and half-integers, i.e., m, I = 0, 4, 1, 1, • • • . 

Let us notice that when there is no purely shear-like potential, i.e., V& = 0, 
then the ^-equation is just nothing else but the eigenequation for the nutation 
^-factor of the stationary states of the spherical top: 

d2 U , . Q d U fm 2 + 2ml cos + I 2 \x \ 

^ + COt ^-( v ^ M* e V U = °- (51) 

The history of this equation traces back to the Reiche-Rademacher theory of 
quantum top |101 1111 127] and to the Wigner theory of irreducible unitary rep- 
resentations of the group SU(2) [T^l [23 [2S] > i-e., roughly speaking, to the one- 
valued and two- valued irreducible unitary representations of the rotation group 
SO(3,K). Then the quantized eigenvalues ee are given by the expression 

2h 2 , 

eej = j'U + 1) 

A 1 

labelled by non-negative half-integer and integer numbers, j = 0, 1/2, 1, 3/2, 
. . ., i.e., j € {0} U (N/2), N denoting the set of naturals. 
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The corresponding eigenfunctions d J m i(Q) were found by Wigner as factors 
in expressions for the matrix elements of unitary irreducible representations of 
SU(2), 

^ m ,($,e,tf) = e im *cP ml (G)e m . 

Here, as mentioned, $, 8, <I> denote the Euler angles parametrization of SU(2). 
Their range is twice larger than the range of Euler angles on the quotient group 
SO(3,R); this is the reason why the half-integer quantum numbers do appear. 

The celebrated functions D 3 m i appear also as stationary states of the quan- 
tized spherical free top. Energy levels are then given by 

^ = ^ + 1), i = Q,±,l,§,-, 

I denoting the main moment of inertia, and of course they are (2j + l) 2 -fold 
degenerate. The labels of basic j-states, m,l, are quantum numbers of projec- 
tions of the angular momentum respectively on the space-fixed and body-fixed 
z-axes: 

h d h d 

Obviously, to, I run over the range + 1, . . . ,j — jumping by one. 

Strictly speaking, in applications concerning the rotational spectra of molecules, 
one has to restrict ourselves to integer values of j, to and I. There are however 
some arguments that perhaps the half integer values might be also acceptable 

Let us also mention that to, I are good quantum numbers also for a more 
general free symmetric top, not necessarily the spherical one. If J, K are two 
main moments of inertia, / doubly degenerate one, then D 3 m i are still basic 
eigenfunctions corresponding to the energy levels 

Eu = —j{j + l) + h 2 (—- —)l 2 . 

They are 2(2j + l)-fold degenerate, namely with respect to the quantum number 
to and to the sign of I. 

One can wonder whether such a symmetric free top in three dimensions, or 
more general three-dimensional top with some external potential, first of all one 
of the shape U(&) (e.g., heavy top), might be useful as a tool for analyzing 
the two-dimensional affinely-rigid body. This is just a question worth to be 
analyzed. 

6 Quantized harmonic and anharmonic vibra- 
tions 

The Schrodinger equations from the previous section may be solved only when 
the explicit form of potential energy is specified. It is clear that simple solutions 
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in terms of known special functions may be expected only when the potential has 
some particular geometric interpretation. For example, this is the case when the 
corresponding classical problem is degenerate and has some hidden symmetries. 

First let us consider the model of the harmonic oscillator potential ([221 . 
Applying the Sommerfeld polynomial method we obtain the energy levels E = 
E a + Ep as follows: 

E=^haj(4n + A+\m-l\ + \m + l\), (52) 

where 

E a = ^-(4n a + 2+\m-l\) 7 Ep = ^ (Anp + 2 + \m + l\) , (53) 

and w = y/C/[i, n = n a + np, n = 0, 1, . . . ,m,l = Q, ±1, .... We may write: 
(i) if |m| > \l\, then m 2 > I 2 and 

E = hiu (2n + 2 ± m) , 

(m) if |m| < |Z|, then m 2 < I 2 and 

E = haj (2n + 2 ± I) , 

(iii) if \m\ = \l\, then m 2 = I 2 and 

E = Hu (2n + 2 ± m) = fiw (2n + 2 ± /) . 

After some calculations we obtain the deformative wave functions f a (a) and 
fp(/3) in the form: 

f a (a) = a a ni + %e- J 5 a2 F 2 (-n Q ; 1 + a; na 2 ) , (54) 

fp{p) = /3-<J + le~%P 2 F 2 (-np; 1 + 7 ; K /3 2 ) , (55) 

where a — h\m—l\, n = y/Cu/h 2 , 7 = ||m + Z|. 

The constant term 4 occurying in the rigorous quantum formula (|52p and 
absent in the quasiclassical one (|24p was in principle expected. This resembles 
the difference between Schrodinger and Bohr-Sommerfcld-quantized harmonic 
oscillators. This is an essentially quantum effect. 

In the classical part we mentioned that the harmonic oscillator model, in 
spite of its academic character, may have some practical utility, and besides, 
it suggests some reasonable anharmonic corrections well suited to certain of its 
degeneracy properties. The mentioned corrections reduce degeneracy in some 
characteristic way and at the same time the model becomes more realistic. On 
the classical and quasiclassical level we discussed the potential (|3"3"|) , i.e., 

Via.0) 
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The model may be rigorously solved on the quantum level and one obtains the 
following formula for the energy levels: 

E = X -hw {*n + 4 + \m + l\ + ^(m - if + . (56) 

The energy in ((55)) depends on an integer combination of the quantum numbers, 
i.e., n = n a + np. The wave functions are as follows: 

f a (a) = a*K-^e-^F 2 (-n a ; 1 + X ; na 2 ) , (57) 

f fi {fi) = (3~< K i + le-^ 2 F 2 (-rip; 1 + 7 ; nf3 2 ) , (58) 

where 



It is seen that the formula for the energy levels is structurally " almost" identical 
with the quasiclassical one (|35p . i.e., 

u ( . , . „ F, 16C/A 



£ = -Tie; [ 4n+ \m + l\ + \j(m - Z) 2 + ,., j 

This is rather typical for systems invariant under " large" symmetry groups and 
based on interesting geometric structures. There is a characteristic shift of 
energy levels, corresponding to the "null vibrations" of the harmonic part of 
the system. Just like on the classical and quasiclassical levels, the system is 
twice degenerate and its energy levels are essentially controlled by two effective 
quantum numbers: n a + rip + \m + l\ and \m — l\. 

Using the formulas (|47|) . ([48]) . i.e., (|4"9"j). ([50]) . we can also quantize the model 

HMD, i.e., 

rr, ns C ( 4 \ 2C - 

The expression for the energy levels E is as follows: 



E = \4n + 4 + \m + l\ + \j (m - if + , (59) 

where n = n r + The functions f r (r), /^(i?) have the form: 

/ r (r) = r-5+ £ K5+ e e~*'' J F 1 2 {-n r ; 1 + 2e; «r) , (60) 
/.?W=(cos0 fsin0 i^i ^-n^,l + ^+ 7 + x;l + x;cos 2 , (61) 
where 



1 / 2u 2C/i 
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For many physical reasons it would be interesting to discuss the model 
(|21j) . however, we were not yet successful in solving explicitly the corresponding 
Schrodinger equation. 

Rigorous solutions for two-dimensional problems may be useful in micro- 
scopic physical problems (vibrations of planar molecules such as S&, CqHq) 
and in macroscopic elasticity (cylinders with homogeneously-deformable cross- 
sections). Applications in dynamics of nanotubes seem to be possible. 

The next important thing to be done is a more comprehensive analysis of the 
status of analogy with Euler angles and the related complexification problems. 
This will be done in a subsequent paper. Some introductory analysis is outlined 
below. 



7 Planar affine body versus spatial rigid body 

It was mentioned above about certain interesting links between mechanics of 
isotropic affine body in two dimension and the dynamics of three-dimensional 
rigid body, more precisely, rigid body with imposed dilatations. Only certain 
analytical aspects, useful in calculations, were stressed there. However, the 
problem is geometrically interesting in itself and has to do with certain com- 
plexification procedures on Lie groups used as configuration spaces. We shall 
analyze this problem in more detail in a forthcoming paper; here we mention 
only a few simple analytical relationships. 

Let us remind that the metric tensor underlying kinetic energy of the planar 
isotropic affine body was given by 



ds 2 = Tr I 



= dx 2 + dy 2 + dz 2 + du 2 



(62) 



the corresponding kinetic energy form reads 
1~~ \ dt dt 



T= ^Tr 




(63) 



where fi denotes the scalar inertial moment. 

For certain reasons it is convenient to use some modified parametrization of 
the two-polar decomposition 



<t> = ODR-\ 

where O, R are proper orthogonal and D is diagonal, namely, 



(64) 



O 



COS tt 



- sin 

* 

cos -j 



D = 



D 1 
D 2 



, R 



sin 4 



C0S \ ™"» 
— sin -|- cos -j 
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and Di = cxp(a + fo/2), D 2 = cxp (a — b/2). It is convenient and instructive 
from the point of view of our analogies to write these matrices as: 



O = exp ( $—(72 



R 1 = exp 



1 



D = cxp a,— (To exp &-03 



1 



where <j v [y — 0,1,2,3) are Pauli matrices; more precisely, er a (a = 1,2,3) are 
"true" Pauli matrices, so 



(65) 





" 1 


" 




' 


1 " 




' -i ' 




1 


^0 = 





1 




1 





, ^2 = 


i 


, o" 3 


-1 



The crucial point for our analogies and links is that the matrices 
T a = 7T (7 a,i a — 1)2, 3, 

are generators of the group SU(2), the universal covering of SO(3,l 
standard commutation rules 

[n,T 2 ] = r 3 , [r 2 , t 3 ] = n, [r 3 , ri] = r 2 . 

And similarly, the matrices 

T\ = in, t 2 = t 2 , r 3 = ir 3 

are generators of SL(2,R) with the standard structure constants, 

[n,T 2 ] = f 3 , [f 2 ,T 3 ] = Tl, [T 3 ,Tl] = -T 2 . 

Obviously, the matrix 

TO = ^0 = 



1 
1 



(66) 
with 

(67) 

(68) 

(69) 
(70) 



generates real dilatations. So, the matrices f v generate the group GL(2,R), the 
configuration space of the planar affine body, and t v generate R + SU(2), the 
2 : 1 covering of the configuration space of rigid body with admitted dilatations 
("breathing top"). The rough symbol ]R + SU(2) denotes the manifold of all ma- 
trices obtained as products of special unitary matrices by positive real numbers, 
M + SU(2) := {\u : A e M+, u e SU(2)}. 

In our models of the doubly isotropic planar affine body, with the met- 
ric element (|62[) we were used rather to parametrize the plane of deformation 
invariants (Di,D 2 ) by r = p 2 — (-D1) 2 + (-D2) 2 and the angle 1} such that 
sint? = (TV - L> 2 2 ) / {Di 2 + D 2 2 ) so that the relationships CE1}-C!I]) and those 
following them are satisfied. However, in models with affinely- invariant kinetic 
energies the variables a, b as deformation invariants are more convenient. As 
mentioned, one can show that 



ds 2 = 



dp 2 + X - P 2 (dO 2 



d$ 2 + 2 cos Qd$dV + rf* 2 ) 
2 (dS 2 + d$ 2 + 2 cos 9d$d* + rf* 2 )) 



(71) 
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We easily recognize the term characteristic for the spherical top described in 
terms of the "Euler angles" ($, 6, ^f) and the term corresponding to the evolu- 
tion of the invariant r, a kind of "dilatation" (not in a rigorous sense). Using the 
more geometric variables a, b and the auxiliary, literally dilatational variable 

6 = VDiD 2 = exp (a/2) , (72) 

we express (fTTj) as follows: 

ds 2 = cosh b dS 2 + S sinh b dSdb +-S 2 cosh b db 2 

4 

1/2 

+ -<5 2 cosh 6 d$ 2 H —d^d^ + d^ 2 

4 \ cosh b 

This is an ugly non-diagonal form; the reason is that ds 2 is not affinely- 
invariant, but only isotropic. The "Euler angles" term is readable, because, as 
we saw, (cosho) -1 = cosG. There are no essential geometric arguments against 
modifying (I71[) by some extra term proportional to dp 2 . 

Let us compare these formulas with those for the spherical three-dimen- 
sional rigid body with dilatations. More precisely, we write down the formulas 
on the group R + SU(2), the universal (2 : 1) covering group of R + SO(3,IR) 
(roughly speaking, the spinorial breathing-rigid-body). Again the rough symbol 
M + SO(3,R) denotes the group of all matrices which are products of proper 
rotations (special orthogonal matrices) by positive real numbers, M + SO(3,K) := 
{XA : A e R + , A G SO(3,R)}. Then </> <= R+SU(2) is " Euler-parametrized" as 
follows: 

4> = exp(aTo)exp($T2)exp(0T3)exp(5'T2). (74) 

More precisely, historical term "Euler angles" is used when the following 
convention is used: 

4> = cxp(aro)exp($T3)exp(9ri)exp(5'T3), (75) 

or similarly, (more popular in textbooks), 

(f>" = exp(a-r )exp(<I>T3)exp(6T2)exp( 1 I'T3). (76) 

If (|74p - (1751 arc identified, then, obviously, ($, O, ^) in those formulas denote 
numerically different functions on SU(2). Nevertheless, there is no essential 
difference between them. What matters is that the SU (2)-matrices are factor- 
ized into products of three elements taken from two orthogonal one-parameter 
subgroups. This is only the question how those three one-parameter subgroups 
are called (ordered). The non- historical, apparently exotic convention (|74l) is 
optimally adapted to our programme of exhibiting some links between planar 
affinc body and spatial rigid body. 

Namely, let us take the following metric on K + SU(2), underlying the kinetic 
energy of the spherical breathing top: 

ds 2 = Tr , (77) 
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where the "f - symbol" denotes Hcrmitian conjugation of matrices. Denoting 
again: 

5 = exp (a/2) , A = S 2 = exp(a), (78) 

we obtain: 

ds 2 = dS 2 + i<5 2 (d6 2 + d$ 2 + 2 cos 6d$d* + d* 2 ) , (79) 
i.e., equivalently, 

ds 2 = — (dA 2 + A 2 (de 2 + d$ 2 + 2 cos Gd^d* + d* 2 ) ) , (80) 
4A 

or, 

ds 2 = ie a (da 2 + d9 2 + d$ 2 + 2 cos 9d$d* + d* 2 ) . (81) 

Obviously, the R + -factor in R + SU(2) is a normal divisor and from the purely 
geometrical point of view of two-side invariant metrics on R + SU(2), there are no 
obstacles against modifying ds 2 by adding an arbitrary correction term ds 2 corr = 
c dS 2 , c being a constant. This means that (|79"|) - (l8T1) may be replaced by 

ds 2 = (1 + c)d5 2 + i(5 2 (d6 2 + d$ 2 + 2 cos 6d$d* + d* 2 ) , (82) 

ds 2 = — ((1 + c)dA 2 + A 2 (d6 2 + d$ 2 + 2 cos 6d$d* + d* 2 )) , (83) 

ds 2 = ie a ((1 + c)da 2 + d9 2 + d$ 2 + 2 cos 9d$d* + d* 2 ) . (84) 

Concerning the extra dilatational term in dynamics of the breathing top, 
cf, e.g., [5]. Replacing the real parameter a in (|77l) . (|78p by imaginary one ia, 
one obtains instead (|84l) the following arc element for the two-side invariant 
Riemannian metric on the unitary group U(2): 

ds 2 = - ((1 + c)da 2 + dQ 2 + d$ 2 + 2 cos 6d$d^ + d* 2 ) . (85) 

For some application or just comparison purposes one can admit in (|77|) . (|78|l the 
general complex parameter a. This results in the doubly-invariant Riemannian 
metric on (C /{0})SU(2) =R+U(2). 

This was, so-to-speak, "one side" of injecting geometry and dynamics of the 
"breathing top" into those of planar affine body (or conversely). There is also an- 
other aspect, namely one based on affinely-invariant metric tensors on GL + (2, K.) 
[22]-[24]. Such metric tensors, of non-definite signature (SL(2, M), GL + (2, R) 
are non-compact, SL(2,R) is semisimple, and GL + (2,R) is the direct product 
of R + SL(2, R)) are linear combinations of those given by the arc element 

ds 2 = Tr (n 2 ) = Tr (f7 2 ) (86) 
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and the purely dilatational correction term 

ds 2 corr = Tr (n f = Tr (o\ 2 , (87) 

where the Lie-algebraic Cartan one-forms D,, ft on GL(2,R) are given by the 
usual formulas: 

Q=(#)<j&"\ h = (f>- 1 d<f> = <p- 1 {l(f>. (88) 

Of course, (|86|) is the main, non-degenerate term of signature (+ + H — ). Killing 
tensor on GL(2,R) is degenerate; the singular direction is that of the one- 
dimensional center R + Id2- This Killing case corresponds to the ratio 4 : (—2) 
of coefficients at {SB), (|871) . 

For calculations we need the following parametrization of cf> G GL + (2,R), 
analogous to (|74)) 

<j> = exp(aTo)exp($T2)exp(6r3)exp(^T2) 

= <5exp($r2)exp(6f 3 )exp(4 r T2), (89) 

where, obviously, _ 

5 = exp (a/2) = vT (90) 

Combining (|86j) , (|87l) with appropriate coefficients (that at the main term (|86[) 
must be non- vanishing) , we finally obtain: 

ds 2 = (1 + c)d5 2 + ^S 2 (db 2 - d<$> 2 ~ 2 coshb - df 2 ) 

= — ((1 + c)dA 2 + A 2 (rf6 2 - d$ 2 - 2cosh6 - rf* 2 )) 

= ie a ((l + c)rfa 2 + rf6 2 -d$ 2 -2cosh6rf$d*-rf* 2 ). (91) 

The relationship between these formulas (as matter of fact, one formula 
written in three alternative forms) and (I7T|) . (|73l) . (|79^) - ([ST|) is obvious. Namely, 
the last four terms in any form of (|9"Tj) become the " minus" terms of the spherical 
top, when some complexification procedure is performed, i.e., when we put 
b = iQ, O being real. Then, obviously, the last four terms become the spherical 
top expression with reversed sign, 

- dO 2 - d<$> 2 - 2 cos 9d$d* - d* 2 , (92) 

and no wonder, because SL(2, R) and SU(2) are two different (and is a sense, 
having opposite properties) real forms of the same complex Lie group SL(2, C). 
The over-all minus term of the Killing metric on SU(2) is due to its compact- 
ness. Performing a similar " imaginarization" of a, we obtain just the "minus" 
expression the doubly invariant metric on U(2). This also expresses the 
fact that GL + (2,R), U(2) are two different real forms of GL(2, C). 
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